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We study the magnetic excitation spectra in L-edge resonant inelastic x-ray scattering (RIXS) 
from undoped cuprates. We analyze the second-order dipole allowed process that the strong per- 
turbation works through the intermediate state in which the spin degree of freedom is lost at the 
core-hole site. Within the approximation neglecting the perturbation on the neighboring sites, we 
derive the spin-flip final state in the scattering channel with changing the polarization, which leads 
to the RIXS spectra expressed as the dynamical structure factor of the transverse spin components. 
We assume a spherical form of the spin-conserving final state in the channel without changing the 
polarization, which leads to the RIXS spectra expressed as the 'exchange'-type multi-spin correla- 
tion function. Evaluating numerically the transition amplitudes to these final states on a finite-size 
cluster, we obtain a sizable amount of the transition amplitude to the spin-conserving final state in 
comparison with that to the spin-flip final state. We treat the itinerant magnetic excitations in the 
final state by means of the l/S-expansion method. Evaluating the higher-order correction with 1/S, 
we find that the peak arising from the one-magnon excitation is reduced with its weight, and the 
continuous spectra arising from the three-magnon excitations come out. The interaction between 
two magnons is treated by summing up the ladder diagrams. On the basis of these results, we 
analyze the La-edge RIXS spectra in Sr2Cu02Cl2 in comparison with the experiment. It is shown 
that the three-magnon excitations as well as the two-magnon excitations give rise to the intensity 
in the high energy side of the one-magnon peak, making the spectral shape asymmetric with wide 
width, in good agreement with the experiment. 

PACS numbers: 78.70.Ck, 72.10.Di, 78.20.Bh, 74.72.Cj 



I. INTRODUCTION 



Resonant inelastic x-ray scattering (RIXS) has recently 
attracted much interest as a useful tool to investigate ex- 
cited states in solids. The K- and L-edge resonances have 
been widely used in transition-metal compounds. When 
one tries to interpret the experimental results, some the- 
oretical support becomes inevitable. Since the interme- 
diate states involved in the scattering process at the K- 
and L-edges are different from each other, the theoreti- 
cal analysis should discriminate the difference with great 
care. 

The K-edge resonance is more suitable than the L-edge 
to detect the momentum dependence of the spectra, since 
the wavelength of the LsT-edge x-ray is an order of the lat- 
tice constant. In the K-edge RIXS, the Is-core electron 
is prompted to an empty 4p state by absorbing photon, 
then charge excitations are created in order to screen the 
core-hole potential, and finally the photo-excited Ap elec- 
tron is recombined with the core hole by emitting photon. 
Charge excitations are finally left with energy and mo- 
mentum transferred from photons. ir— The spectra have 
been analyzed by several methodsi^^— Among them, a 
formalism developed by Nomura and Igarashi (NI)i2i"— 
usefully describes the RIXS spectra in terms of the Sd- 
density-density correlation function on the initial state 
It is based on the Keldysh Green functionji^ and is re- 



garded as an adaption of the resonant Raman theory of 
Nozieres and AbrahamsJ^ In the actual application of 
the formalism, the electronic states are described within 
the Hartree-Fock approximation on the multiband tight- 
binding model, and the electron correlation is treated 
within the random phase approximation. The spectra 
are calculated as functions of energy loss and momen- 
tum transfer in good agreement with the experiments for 
undoped cuprates, ^^^^^^ NiO,^^ and LaMnOs.^^ 

In addition to the charge excitations, the magnetic 
excitations have been observed from the Cu LC-edge 
RIXS experiment in La2Cu04j2£!31. They are brought 
about through the modification of the exchange cou- 
pling between the spins of 3d electrons at the core- 
hole site and those at the neighboring sites by the 
core-hole potentiali^"— To investigate the magnetic ex- 
citation spectrum, the NI formula has been adapted 
in our previous worki^ By replacing the channel of 
creating an electron-hole pair by that of creating two 
magnons, the expression is derived that the RIXS spec- 
tra are proportional to the two-magnon correlation func- 
tion. The magnon-magnon interaction in the correla- 
tion function has been treated with the use of the 1/5- 
expansion method [S is the magnitude of spin)j^"— The 
obtained spectra, reflecting a significant influence from 
the magnon-magnon interaction, show characteristic de- 
pendence on the energy loss and the momentum trans- 
fer, which captures features the experimental data have 
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demonstrated ] 1 

Unlike the iC-edge RIXS, the L-edge RIXS directly 
accesses the 3d states, which leads us to a necessity to 
develop an appropriate treatment for the L-edge RIXS. 
The process is illustrated for undoped cuprates in Fig. [T] 
the 2p-core electron is prompted to the empty — 
orbital by absorbing photon, and then an occupied 3rf 
electron combines with the core hole by emitting photon. 
Here x, y, and z axes correspond to the crystal axes a, b, 
and c, respectively. If the 3d state in the photo-emitting 
process is different from the — y^ orbital [Fig. [T] (a)], 
excitations within the 3d states are left in the final state, 
which process is called the 'd-d' transitioujS^ Since such 
excited states are quite localized at the core-hole site, the 
spectra show little momentum dependence. In this paper, 
skipping the study of the d-d transition, we concentrate 
on the study of the L-edge RIXS spectra arising from 
magnetic excitations, which have recently been observed 
in undoped cuprates<^"— 

It is known that the spin- flip excitations [Fig.[TJb)] give 
rise to the spin-wave-like dispersion, when the direction 
of the staggered moment deviates from the z direction. '^^ 
Since the spin angular-momentum is coupled to the or- 
bital angular-momentum through the large spin-orbit in- 
teraction of the 2p-core states, the change by the spin- 
flip may be compensated by the change of the polariza- 
tion of photon. In addition, it is expected that the spin- 
conserving excitations [Fig. [Ijc)] could be brought about 
through the intermediate state, since the strong pertur- 
bation is working by losing the spin degree of freedom at 
the core-hole site in the intermediate state. 

For investigating the magnetic excitations, the RIXS 
spectra have been analyzed by assuming an extreme con- 
dition that the core-hole life-time is so short that the in- 
termediate state can not have enough time to relax. 
Only the one-magnon excitation is brought about with- 
out any two-magnon excitations. This is called as the 
ultra-short core-hole life-time (UCL) approximation,-- 
which reminds us of the fast collision approximation in 
resonant elastic x-ray scattering i^i^^ In reality, the spec- 
tra observed as a function of energy loss exhibit a sys- 
tematic change with changing momentum transfer;'^- the 
peak position moves according to the spin-wave disper- 
sion curve, while their shapes exhibit structures indica- 
tive of two- or three- magnon excitations. Accordingly the 
assumed condition is obviously unsatisfied in undoped 
cuprates. The purpose of this paper is to develop a com- 
prehensive theory describing the spectra beyond the UCL 
approximation. 

The situation that the spin degree of freedom is lost at 
the core-hole site in the intermediate state is similar to 
the situation of the non-magnetic impurity, which prob- 
lem has been extensively studied by the linear-spin- wave 
(LSW) methodi^i^ The method developed there, how- 
ever, is not applicable to the present problem, since the 
calculation of the RIXS spectra require not alone the in- 
termediate state but also the initial and final states in 
evaluating overlaps between them. It may not be logi- 
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FIG. 1: La-edge RIXS process for undoped cuprates. Black 
and white circles represent electrons and holes, respectively. 
Vertical arrows indicate the spin direction. Wavy lines indi- 
cate the incident and scattered photon. The t2g orbitals are 
omitted. 

cally appropriate to perform a perturbation calculation 
with the terms involving the lost spin degree of freedom, 
although such attempts have been mentioned and tried4i 
Along this consideration, the perturbative expansion like 
the NI formula may not be sufficient to treat the inter- 
mediate state. 

In this paper, we analyze the second-order dipole al- 
lowed process that the strong perturbation works through 
the intermediate state. Within the approximation that 
the perturbation is not extending to neighboring sites, 
we derive the spin-flip final state expressed as OLf± x 
f^i-L • S0I5) in the scattering channel with changing the 
polarization, where cti(/)j^ is the polarization vector of in- 
cident (scattered) photon projected onto the x-y plane, 
and So is the spin operator vector at the core-hole site. 
This form is inferred from the formula of the elastic 
scattering.— 1^ This final state leads to the RIXS spectra 
expressed as the dynamical structure factor of the trans- 
verse spin component. In the scattering channel with 
changing the polarization, we assume a spherical form of 
the spin-conserving final state, which leads to the RIXS 
spectra expressed as the 'exchange'-type multi-spin corre- 
lation function. In the evaluation of the transition ampli- 
tude to these final states, we introduce a finite-size cluster 
and use the wave- functions numerically obtained by diag- 
onalizing the Hamiltonian. We obtain a sizable amount 
of the transition amplitude to the spin-conserving final 
state in comparison with that to the spin-flip final state. 

We express the spin-flip and spin-conserving final 
states in terms of magnon creations and annihilations 
with the systematic use of the 1/S expansion. The spin- 
flip operators are expanded in terms of one-magnon 
creation and annihilation operators within the LSW the- 
ory, and of three-magnon creation and annihilation op- 
erators in the second order of 1/S ^"—1^1^ Hence the 
spin-flip final state gives rise to the RIXS spectra con- 
sist of two ingredients. One is the (5-function peak with 
its weight reduced from the LSW value and the other is 
the broad three-magnon continuum as a function of en- 
ergy loss. The spin-conserving operators are expanded 
in terms of two-magnon creation and annihilation op- 
erators with neglecting the higher-order terms such as 
four-magnon creation and annihilation operators. Since 
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two magnons are created near the core-hole site, we take 
account of their interaction by summing up the ladder di- 
agrams, as was done in the magnetic excitation spectra 
in the if-edge RIXS^ 

On the basis of these results, we analyze the Cu 
La-edge spectra of Sr2CuCl202, which show strong de- 
pendence on the polarization in accordance with the 
experiment)^ We find that a considerable amount of 
satellite intensity comes from the three-magnon excita- 
tions in the spin-flip final state, which leads to the asym- 
metric shape as a function of energy loss in the tt polar- 
ization. We also obtain the substantial contribution from 
the two-magnon excitations in the a polarization. These 
results agree well with the experiment. 

The present paper is organized as follows. In Sec. 
we describe the Hamiltonian responsible for magnetic ex- 
citations and transition-matrix elements relevant to the 
^2,3 edges in cuprate compounds. In Sec. IIIH we dis- 
cuss the process giving rise to the magnetic excitations 
through the intermediate state. The RIXS spectra are 
expressed in terms of spin-correlation functions. In Sec. 
IIV| spin-correlation functions are calculated by means of 
the 1/5* expansion. In Sec. El the calculated Ls-edge 
spectra of Cu are compared with those for Sr2Cu02Cl2. 
Section |VT] is devoted to the concluding remarks. In Ap- 
pendix, a brief summary of the l/S'-expansion method is 
given. 



II. FORMULATION OF RIXS SPECTRA AT 
THE L2,3 EDGE 

A. Hamiltonian 

At the half-filling in cuprate compounds, each Cu atom 
has one hole in the — orbital, where the x and y 
axes are defined along the Cu-0 bonds and the z along 
the crystal c axis. Adopting the hole picture, we assume a 
single band Hubbard model on a two-dimensional square 
lattice for 3c? electrons: 

H = tY. {did,, + H.c.) + uY,d\^d\^d,^d,^, (2.1) 

where di, (c?Jg.)represents the annihilation (creation) op- 
erator of the hole with spin a at site i. The sum over {i, j) 
extends over distinct pairs of nearest neighbors. Hopping 
integral and on-site Coulomb interaction are denoted as 
t and U, respectively. We have neglected the small inter- 
layer coupling. This Hubbard model may be mapped 
from a more precise "d-p" model for cuprate compounds. 

For low-energy spin excitations, the Hubbard Hamilto- 
nian is further mapped onto the two-dimensional Heisen- 
berg Hamiltonian with the exchange coupling constant 
J = 4t'^/U, 

^mag — J ^ ^ Si • Sj. (2-2) 



This Hamiltonian leads to an antiferromagnetic order. 
It is known that Cu spins align on the Cu02 plane at 
the angle of n/A to the Cu-0 bonds.— Therefore, it is 
convenient to define the spin coordinate frame x', y' , z' 
axes by rotating the crystal-fixed coordinate frame a, b, 
c axes through the Euler angles a = 7r/4, (3 — 7r/2, and 
7 = 0. The t and J, states in Eq. (|2.ip are interpreted as 
the eigenstates with respect to . Hereafter, we shall 
describe the spin state of the 3d electron in the x'y' z' 
coordinate. 



B. E\ transition at the 1/2,3 edge 
The Hamiltonian of photon may be written as 

where Cq^ ^A^ii) stands for the annihilation (creation) op- 
erator of the photon with momentum q, energy Wq, and 
polarization direction [i {— x^y, and z). In the electric 
dipole {El) transition, a 2p-core electron is excited to 
the 3c? states at the transition-metal £2,3 edge. The 2p 
states are characterized by the total angular momentum 
j = 3/2 and 1/2 due to the strong spin-orbit interaction. 
The eigenstates with j — 3/2 may be expressed as t)j 

yi73|<^i i) + v^l^o t), ^273100 i) + yT73|^-i t), 
\(f)-i i), for m = 3/2, 1/2, —1/2, —3/2, respectively, and 
those with j = 1/2 may be expressed as — ■y/2/3|(/)i I 

) + yT73|0o t), -\/T73|0o i) + v/273|</'-i t), for m = 1/2 
and —1/2, respectively, where m represents the magnetic 
quantum number. The orbitals 01, 0o, and 0_i have an- 
gular dependence Yn, Yio and Yi_i, respectively. In the 
above expression, the coordinate frame for spin is defined 
the same as the orbitals, which is different from that for 
the spins of the 3c? states, that is, f and are associated 
with the direction of the crystal c axis. Therefore, the 
interaction between photon and electron at site i may be 
described as 

Hint^wJ2^= E D^{jm,a)hl^c^^d,,e'^''^+R.c., 

q,/j. V ^ ■i,m,a 

(2.4) 

where h^^^ stands for the creation operator of the 2p hole 
with the angular momentum jm. The w is a constant 
proportional to R3d{r)R2p{r)dr, with R3d{'r) and 

R2p{f) being the radial wave- functions for the 3c? and 
2p states of Cu atom. The D'^{j'm,,a) describes the de- 
pendence on the core-hole state and 3c? spin, which is 
calculated by taking care of the difference in the defini- 
tion of the spin axes between the 3c? and 2p states. Table 
H] shows the calculated values of D'^{j'm, a). Note that no 
El transition takes place for the polarization parallel to 
the z axis. 
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TABLE I: Numerical values of D''{jm,a) 



Here, 



g-v- i '"-cos^ and v = e*^^""" sin^, where a, fi, and 7 
are the Euler angles transforming the crystal-fixed a, b, and 
c axes to the spin axes x' , y' , and z' axes. 



J 



(7 =4- 



3/2 



1/2 



3/2 
1/2 
-1/2 
-3/2 
1/2 
-1/2 



C. Formulation of RIXS spectra 

Following Nozieres and AbrahamSfi^ we use the 
Keldysh-Schwinger formalism^'' to investigate the RIXS 
spectra. First we prepare the initial state 



1$,; 



(2.5) 



where \g) represent the ground state of the matter with 
energy Eg. The incident photon has momentum q^, en- 
ergy Wi, and polarization direction a^. Then we calculate 
the probability of finding a photon with momentum q^, 
energy lo / , and polarization direction a/ at time to using 
the following formula, 

(2.6) 

Here U{t,t') is the S matrix defined by 

U{t, -cx)) = Texp |-i J ffi„t(t')dt'| , 



(2.7) 



where 7?i„t(t) = exp[i{H + Hph)t]Hiy,texp[~i{H + HpY,)t] 
with H and Hp\^ being the Hamiltonian of matter and 
photon, and T represents the time-ordering operator. 
The S-matrix is expanded up to the second order with 

Hint ■ 



U{t, -00) = 1 + 

+ (-0 



Hint{t')dt' 



t pt' 

2 / dt' 



dt"Hintit')Hintit"). 

(2.8) 



By substituting this into Eq. (|2.6p . we obtain the 
transition probability per unit time with to — > 00, 



W{qfaf; qiai) in the following form: 

W{qfaf\q.iai) 

/O nOO nU 

dt du' di'^^($,|i7i„t(t')l^'> 
-00 J-00 J-00 J, 
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3/2 










1/2 










-1/2 










-3/2 




Vio 




1/2 


1/2 






X 




-1/2 










f 



int\n) 



E„ + LLJi — En 



(ri|i?i„t|$0 



X 5{Eg +uji~ Ef, - LOf), 



(2.9) 



with qi = (q.j,Wi), g/ = (q/,w/), and |$/.) = Cg^a/I/')- 
This is nothing but conventional expression of the second- 
order dipole allowed process, on which our following anal- 
ysis is based. The |/') represents the eigenstate of the 
Hamiltonian of the matter with the eigenvalue Ef , while 
|n) represents the eigenstate with eigenvalue £"„ in the 
presence of core hole. 

In our previous paper having discussed RIXS at the Cu 
K edge^ the E\ transition could not directly change the 
3d states but simply create the core-hole potential which 
attracts electrons in the 3c? states. In such a situation, 
by treating the core-hole potential as a perturbation, we 
could conveniently introduce the Green functions and di- 
agrams on the basis of the Keldysh-Schwinger formalism 
to calculate W{(qfaf \ qtai). To the present case, however, 
such diagrammatic procedure is difficult to apply, since 
the El transition modifies the 3d states from the 
configuration to the 3(i^*^-configuration, which change is 
hard to express by the diagrams for magnon excitations. 



III. MAGNON EXCITATIONS AROUND THE 
CORE-HOLE SITE 

A. General consideration 

Assuming that the core hole is created at the origin 
in the intermediate state, we analyze the spin system 
around the origin, as shown in Fig. [2] We write the 
ground state \g) of H^ag as 



\9) = \ml) + \m'o), 



(3.1) 



where | 'f) and | l) represent the spin states at the origin, 
and I'^q) and |?/'o) are constructed by the bases of the 
rest of spins. The expectation value of at the origin 
may be expressed as 



(3.2) 



Since the spin degree of freedom is lost at the origin 
just after the El transition takes place, the wave function 
becomes 



Hintlg) 
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J2 D"^ijm,a)\^^] 



\jm), (3.3) 



5 



o 



o- 



-0- 



-o 



■o- 



■o 



7 4 



-o 



o- 



■o- 



■o 



o 



FIG. 2: Spins near the core-hole. Each site is labeled by 
number and the core-hole site is assigned number 0. 



where \jm) represents the core hole state. Let H' be 
the Hamiltonian in the intermediate state. It consists 
of the exchange interaction terms excluding those with 
the core-hole site. The Hilbert space representing H' is 
different from the initial state. This situation is the same 
as the problem of the non-magnetic impurityj^"^ Note 
that the states and |-0g) are not the eigenstates of 
H' . They are expanded in terms of normalized eigenstate 
|(/)^)'s of H' with eigenvalue ej^. We therefore have 



oc 



with 



n ^ 

J2W)\(b,)R{e'^){^,\ro'), 



R{e') = 



(3.4) 



(3.5) 



where Cg represents the ground state energy of 7?mag- 
The ecore denotes the energy required to create a core 
hole in the state \jm) and the 3d^°-configuration. The 
r stands for the life-time broadening width of the core 
hole; r ~ 0.3 eV at the L3 edge. Since H' is different 
from iJinag only around the core-hole site, \4>n) could have 
sufficient overlap with ) only when — Cg varies in 
the range of several J's (J ^ 0.1 eV). Note that the spin 
degree of freedom at the origin is recovered in the final 
state and the wave function is given by |cr)|(/)^). The first 



factor in Eq. (|3.4[) has relations, 

Y^D^f{jm,aYD'-^{3m,a) = p(°)(j; a/, a0(3.6) 

:i^"^(jm,a)*I?"'(jm,-a) = Pi^^Cj; «/, a,)(3.7) 



Here, for a ~a denotes 4, and vice versa. Table HIl 
shows pjP^ and Pi^"* for ai and af along the x, ?/, and z 
axes. The extension to the cases of general and a/ di- 
rections is obvious. The Pi^"* and P^^' correspond to the 
spin-conserving and the spin-flip processes, respectively. 
Note that Pi^^(l/2; a/, a,) Pi^^(3/2; a/, a,) = 0. As 
will become clear later, spin-flip excitation spectra are 
proportional to jP'^-'^^p. Therefore, if the processes for 
j = 3/2 and j — 1/2 are not separated, no spin- flip exci- 
tation comes out. 



B. Scattering channel with changing the 
polarization 

We analyze the scattering channel that the polarization 
changes; let ai and a/ be along y and x axes, respectively. 
We assume that the staggered magnetic moment is par- 
allel to the z'-axis, where the coordinate frame of the 
x\y' , z' axes is specified by the Euler angles a, /3, and 7 
from the coordinate frame of the x, z axes. We have 
the spin-conserving term coming from P^"'' and the spin- 



flipping term coming from Po- 
term is given by 



(1) 



1 



Eg — En 



The spin-conserving 
(n|i/int|g) 



oc 



15 



cos/3 |t)5]|0,)P(e;)(0,|V4) 

;>Ei'^'')^K)('^''i^o)|- 



(3.8) 



Since this state has an overlap to the ground state, we 
have the elastic amplitude Adas, 



Adas (X (^--j cos/3 |(V'J|^|0,)P(e;)(0„|4) 

-(^o^l5^|0,)P(e;)(0,|^J[)l. (3.9) 



Introducing the quantity 
1 



{^S\J2\cj,^)Rie'^){q^,\ro), (3.10) 
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TABLE II: P^°^ 




) and Py 


(j;af,ai) where upper and lower signs 
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with a =t and I, we define fQ^\u!i) and A{uji) by 

/W(..0 = /^'^M±^A(a.,), (3.11) 

where plus and minus signs in the second term cor- 
respond to a =t I, respectively. In the far-off- 
resonance condition that Wi <^ Ccoro or w.^ ^ ecorc, and in 
the large limit of F, which is called as the UCL condition, 
that r ^ \uji — Ecorc I, |ej, — Egl, we could factor out i?(e^) 
from the summation over 77 in Eq. (j3.10l) . Thereby, using 
Si; \'t>v){4'v\ — 1' immediately obtain A(wi) = 0. 

By inserting Eq. p. lip in Eq. (13. 9p . and by using Eq. 
(13.21). we have 



^cias (X (^-— j COS/3 |2/^^^(c^0(5'o ) + 2^1^^ 

(3.12) 

This result may be compared with the conventional ex- 
pression of the elastic magnetic scattering amplitude by 
Hannon et al.,— 



(T^^^a/ X ai ■ m, 



(3.13) 



where m is the staggered magnetic moment vector, and 
cr'-^'' is a certain numerical constant. In order that 
Eq. (|5T^ is consistent with Eq. (|XTO)) . A(a;j) could be, if 

it existed, expanded a.s a J2s i'^s )+^J25'{^S')~^' ' ' '^i^^ 
6 and 5' denoting the nearest and next nearest neigh- 
bor sites to the core-hole site (a and b have to go 
to zero in the far-off-resonance condition as well as in 
the UCL condition), since A(wi) should be incorporated 
into the renormalization of in Eq. (j3.13p . The pres- 
ence of A(cLii) may suggest that the disturbance through 
the intermediate state reaches to the neighboring lat- 
tice sites. Actually, Eq. p. 81) has finite overlaps to 
Sg \g) and 5|, \g) , which are not orthogonal to both Sq \g) 
and 5*0 By neglecting such possibility, we project 
Eq. p.8p onto \g) and Sq \g). Since \g) is not orthog- 
onal to \g) and not normalized, we need to introduce 
overlap matrix p between \g) and 5'q \g). Using the in- 
verse of yO, we have 



Hi, 



En — En 



{n\Hint\g) 



= (-^) cos/3 |iA(a;,)|5) + 2f^'\u;.)S^' \g){^..U) 



By neglecting A(aji), the final state in this channel is 
approximated as 



cos/32/«(..,)5o^'|g) 
^) 2f^^\uj,)oLf^ X - Souls), (3.15) 



where Sq[| stands for the component of Sq parallel to 
the direction of the staggered magnetic moment. The 
OLi±^ and a f± represent the polarization vectors projected 
onto the x-y plane. Therefore afj_x is always paral- 
lel to the z axis. Note that the inelastic terms are some- 
times inferred from Eq. (I3.13P with simply replacing m 
by the spin operator So at site 0.— Equation p.lSp is, 
however, different from such a term, since Eq. p.isp is 
restricted within the spin-conserving process, and disap- 
pears for /3 = 7r/2. 

The final states responsible to inelastic scattering 
mainly come from the spin-flip terms. 



1 



int\n) 



Eg — En 



(nli/i„tl$.) 



(X pf)|;)5]|0,)i?(e;)(0,|Vo^) 
n 

i)^|</.,)i?(e;)(0,|Vo^) 



cx ( -)sin/3<;e-*^ 



p*7 



t)^|<^,)i?(e;)(0„|Vi) 



(3.16) 



We project this state onto S^lg) with Sq = Sq ± iSQ . 
These states are orthogonal to each other and to \g) 
but not normalized, that is, {g\SQ SQ\g) — (V'olV'o) 
and {g\SQ SQ\g) = (V'olV'o)- Thereby we could express 
Eq. ([XTC]) as 



sml3{f^'\u,){e-^^S-+e^^S^)\g) 
+ ^A{uj,){e-^^SQ-e^^S+)\g)}. (3.17) 
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The first term may be rewritten as 



with 



sin 132 f^'\uj,){cos^SS' -sin jSi)\g) 
2f^^\uj,)af^ X a,^ ■ So±|.g), (3.18) 



where So± represents the component perpendicular to 
the direction of the staggered magnetic moment. As 
regards the second term of Eq. (I3.17p . the inclusion of 
A(LiJi) would require adding the states S^\g) and Sf;\g) 
as the states to be projected. Such an analysis would be 
rather complicated and will not be attempted in this pa- 
per. With disregarding the term of A(a;i), we finally have 
the expression of the final state by combining Eq. p.lSp 
with Eq. dXTT]) : 



E 



1 



Eg — En 



•Sols) 



(3.19) 



C. Scattering channel without changing the 
polarization 

In this scattering channel, only the spin-conserving ex- 
citations are brought about through the diagonal com- 
ponents of For cx.i — OLf — (1,0,0) and for 
OLi ~ ocf — (0, 1,0), we have 



(^i + E„ — En 



{n\Hint\g) 



;)5]|0,)i?(e;)(0,|Vi) 



(3.20) 



We consider the spherical form of the final state, X • 
S0I5) with X = Sg. The deviation from the spherical 
form would be considered as a next-step approximation, 
since we are neglecting A(aJi) which suggests that Ss\g) 
are to be included as possible excited states. Since X • 
S0I5) is not orthogonal to \g), we introduce the overlap 
matrix p defined by {p)ij = {ipi\'ipj) with = \g) and 
\7p2) = X • S0I5). Then, using the inverse of p, we project 
the final state onto these states, resulting that 



1 



cx /f)(w.)l5) + /r(^.)X-So|5) 
where /m (i^i)'s are given by 



Eg — E„ 

(2) I 



{n\Hint\g) 



(3.21) 



(3.22) 



Q2\^^) - l{Y.sgn{a){ro\X''Y,\cj^,)Rie'^){^,\ro) 

+(V'J|x-^l0„)i?(6;) (0,1^0^) 

+ (^oV+El'^'')^«)('^''l'^^)}' (3.24) 



where sgn(CT) = 1 for a =t and —1 for a =|. 

To understand the possible spin-flip at the origin in 
Eq. p.2ip . we take up the second term of Eq. p.24p . 
and examine the associated generating process. Con- 
sider that the down spin at the core-hole site is anni- 
bilated by absorbing photon. The wave function for 
surrounding spins is expressed by l-^o): which satisfies 

{ti}i\{Y.rSI)\il^i) = 1/2 with i running aU the lattice 
sites except the origin. At the end of time evolution 
in the intermediate state, this wave function is modified, 
and could have a finite overlap with the state X+jV'o), 

since (V'ol^"(E» 'S'f')^+|V'o) = 1/2. Then, the 3d hole 
with down spin is created at the core-hole site by anni- 
hilating the core hole with emitting photon, leading to 
I DA'+IV'J). Since SQ\g) = \ DIV'q), this state is nothing 
but X+S^\g). 

In the far-off-resonance condition and in the UCL ap- 
proximation, i?(e^) could be factored out in Eqs. p.23p 
and (|3.24p . It leads to Qm{^i) — ^(eo)(p)rn,i and as a 
consequence. 



'{P)r 



OJ'Jm,!- 



(3.25) 

Therefore, in both cases, the spin-conserving excitations 
could not be generated. 

One may formally write the intermediate state Hamil- 
tonian H' by eliminating the bonds which connect the 
spin at the origin to spins at neighboring sites from the 
initial state Hamiltonian H-axa„: 



H' — H„ 



(3.26) 



The Hilbert space representing Eq. p.26p formally con- 
tains the spin-degrees of freedom at the core-hole site, 
which should be completely decoupled from the outer 
spins in the final solution, as known from the non- 
magnetic impurity problem.^>i2, Finite-order perturba- 
tion with V could not satisfy this criteria. In this context, 
although a first-order perturbation with V has been at- 
tempted to include two-magnon excitations by extending 
the UCL approximation,**^ it may be logically inappro- 
priate. 
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D. Cluster model 

It is not easy to evaluate accurately fQ^\uJi) and 

(2) 

/m (a;i)'s. In this paper, we use a cluster consisting of 
a central spin and of 8 neighboring spins, as shown in 
Fig. [2] The outer spins labeled as Si ~ Ss are assumed 
under the staggered field from spins outside the cluster. 
For the central spin Sq belonging to the A sublattice, the 
initial-state Hamiltonian is given by 

4 8 

i/mag = ^ E ■ + 2 ^ - 2 2 E ^1' (3.27) 

(i.J) 1=1 3=5 

which is represented by the matrix with 512x512 dimen- 
sions. The intermediate-state Hamiltonian is given just 
by eliminating So from Eq. p.27|) . which is represented 
by the matrix with 256 x 256 dimensions. Diagonaliz- 
ing numerically the Hamiltonian matrices, we obtain the 
eigenstates of iJmag and H' . We obtain, for example, 
(5^) = 0.394 from Eq. (|X2|) . This value is an overesti- 
mate of the sublattice magnetization in comparison with 
the l/S'-expansion value in second order 0.307,'^*' prob- 
ably due to the finite-size effect. Using the eigenstates 
of TJmag and H' , we evaluate /q^^ from Eqs. p.lOj) and 

(|3ll|) . and f^^ from Eqs. and ((STMt . Since the 

cluster size is strongly subject to the boundary, the values 
thus obtained may be considered as a semi-quantitative 
estimate. 

Using the notation in the preceding subsection, we ex- 
press the L2,3-absorption coefficient Aj(uji), 



{uJi + eg - ecorc - e'j^Y + 

Substituting the eigenstates of the cluster model into 
Eq. (|3.28p . we calculate Aj{uji). Figure [3] shows the cal- 
culated Aj{uji) as a function of photon energy. The ori- 
gin of photon energy is uJi = ecore- The dimensionless 
life-time broadening width of the core hole is chosen as 
r/(2J) = 1.2. The calculated curve is found very close 
to the Lorentzian shape. By comparison with the exper- 
imental curve for the Cu Ls-edge in Sr2Cu02Cl2, F is 
estimated ^ 2. 4 J with J = 130 meV. 



E. Spin correlation function 

The core hole could be excited at all the Cu sites in 
the x-ray scattering event. For the contribution from the 
core- hole site r^, we need to multiply weight exp(zq • rg) 
(q = Qi— q/) to the amplitudes discussed in the preceding 
subsections. Collecting the contributions from each core- 
hole site in Eq. (|2.9I) . we finally obtain the expression of 




-10 -5 5 10 

(n/(JSz) 



FIG. 3: Absorption coefficient A(u}i) as a function of photon 
energy uji. r/(2J) = 1.2. The origin of energy is set to cor- 
respond to UJi — ecoro- The dotted line represents the XAS 
experimental data at the Cu L3 edge in Sr2Cu02Cl2,— where 
J is assumed to be 130 meV, and the experimental curve is 
shifted such that the peak position coincides with the calcu- 
lation. 

the RIXS spectra: 

+ |p(")(j;a/,a.)pr(2)(c.,;q,w)},(3.29) 

where the spin-flip and the spin-conserving correlation 
functions are defined by 

/oo 
(Z(i)t(^.q,i)^(i)(^.q,0))e*'^*dt, 
-00 

(3.30) 

/OO 
(Z(2)t(,,^;q,t)Z(2)(^^;q,0))e*'-*di, 

(3.31) 

where the angular bracket denotes the expectation 
value in the ground state. Operators .Z(^^(wi;q) and 
Z^'^\ijji; q) are given by 

Z(i)(c.,;q) = 2fl,'\u,){st [-ci) + St' i-ci)] (3.32) 
Z(2)(c.,;q) = f^^\oj,)[{X^'S^')a{-<l) + {X^'S^'),{-Ci) 

+ (X-'S^')a{-ci) -f (X-'5-')b(-q) 

+ (X«'5'«')a(-q) + (X'^'5^')b(-q)}. (3.33) 

The Fourier transforms of and ^Sq are 

defined separately for the A and B sublattices and dis- 
criminated by subscripts a and 6, respectively. For ex- 
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(a) 



(c) 



D 
X- 



D 




y(^) (wi; q, cj) is proportional to the "exchange" -type 
spin-spin correlation function discussed in the magnetic 
excitations in the JsT-edge RIXS.— 





(d) 



FIG. 4: Diagrams for the time-ordered Green function 
y'^' (a;^; q, a;) required to evaluate the correction up to the 
1/(25)^. Solid lines represent the unperturbed Green func- 
tions GL°2(k,a;) or G^''^(k,a;), on which the arrows and the 
type {a or (3) are omitted. The crosses in diagram (a) rep- 
resent the spin reduction factor. Solid circles represent the 
three-magnon terms given by Eq. (|4.ip and their Hermitian 
conjugates. 



IV. SPIN CORRELATION FUNCTION WITHIN 
THE 1/5 EXPANSION 

Spin-flip and spin-conserving excitations derived in the 
preceding section could freely propagate in the crystal 
in the final state because of the absence of core hole. 
Therefore, the cluster model with small size would not 
work well in the final state. We exploit the 1 / S'-expansion 
method in the study of the final state. Notations and 
several relations required in the present study are briefly 
summarized in Appendix. 



ample, by multiplying weight e"i ""', we have 

(X-'^-\(-q) = ^i5^^f;5^f'e'''-'-% (3.34) 

ieA s 



jeB ^ s 



Note that F^^^ (cj.;; q, w) is proportional to the dy- 
namical structure factor of the transverse spin, while 



A. Spin-flip excitation spectra 

We express the spin-flip operator Z^^\uji;c[) in terms 
of magnon operators by using Eqs. (jAl|) - (jAlO[) . Since 
Z(^^(wi;q) is proportional to (~q) + (— q), the 
derivation is parallel to those of Eqs. (6.9) and (6.10) 
in Ref. 45 in the study of the dynamical structure factor. 
In the following, we simply write down the result. 

For q being inside the first magnetic Brillouin zone 
(MBZ), we have 



2 

1 2 

25 7V 



E'5G(q 



234 



4)-V244 



aLga^ + /3-2a3a4 + sgn(7G )(4/3l3/3l4 + a_2/33/34)}] ,(4.1) 



where 

^q234 = -a;2 + sgn(7G)a;qX3a;4, (4.2) 

(2) 

^q234 = X3a;4 - sgn(7G)a;qa;2. (4.3) 

The Kronecker delta i5g (q + 2 — 3 — 4) indicates the con- 
servation of momenta within a reciprocal lattice vector 
G and sgn(7G) denotes the sign of 7g. is defined by 
Eq. (fATT|) . The spin reduction factor D, which is related 
to the zero-point reduction of spin, is given by 

A5_1A5(1 + 3A5) 

25 4 (25)2 ' ^ > 



with 

k 

For the square lattice, AS = 0.1966. 

Introducing the time-ordered Green 

function defined by ^(^^(wi; q, w) = 
-i/(T(2'(i)t(t^..q^t)^(i)(^..q^0))e'"*dt, we expand 
the Green function up to the second order with 
1/5, which diagrams are shown in Fig. |4l We ob- 
tain the correlation function from the time-ordered 
Green function by the fluctuation-dissipation theorem, 
F^-'^^ (wi; q, w) = —2lmY'' ^\u}i;q,u}). Here, ImA means 
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(0,0) 



Present work (LSW) 

--- UCL 




(it,0) 



(lt/2,7l/2) 



(0,0) 



FIG. 5: (Color online) One-magnon intensity 

yi{uji; q) /[N{2S)] within the LSW contribution, as a 
function of q along symmetry lines, in comparison with the 
results of the UCL approximation (red broken line). The uii 
is set to give rise to the peak in the absorption spectra. 



imaginary part of quantity X. The correlation function 
is found to consist of the (5-function peak of the one- 
magnon excitation and the continuum of three-magnon 
excitations, that is. 



y(i)(wj;q,cj) = yi(cj,;q)(27r)(5(w- eq) +2/3(a;j;q,w). 

(4.6) 

Within the second-order of 1/S, the one-magnon spectral 
weight yi (q) may be expressed as 



2/i(w,;q) = Ni2S)\f['\LU,)m^ + f^'\uj,)e^ 

1 + + (^yi.2(q) !^ .(4-7) 



where 



2/1, i(q) 
2/1,2 (q) 



with 



1 

'2' 



(4.8) 



-iA5(l + A5) + — ES(q,eq) 



2 \ 

) 2^ ^<^q*-pV^q+p-p' 

PP' 



2 \ 

I \ " 9/;2 n2 jfl ij2 

]\J J ^*^q<p*p'<^k+p-p' 

PP' 



-IB 



(4) 

q,p,p',[q+p-p'] 



■^q+p-p 



^ ' q,p,p ,[q+p-p ] ' 

)^ (eq + ep-f ep'-f eq+p_p/)2 



B 



(4) 



(1) 



q,p,p',[q-l-p-q'] q,p,p',[q+p-p 

(fq ~ Cp ~ £p' — £q+p-p') 



sgn(7G)B 



(6) 



r(2) 



q,p,p',[q+p-p'] q,p,p',[q-l-p-p'] 



q,p,p'[q+p-p'] q,p,p',[q+p-p'] 



N 



(2) 

q,p,p'[q+p-p'] 



M 



(1) 

q,p,p',[q+p 



sgn(7G)M^^^ ^ „ 

6 1^' q,p,p',[q+p-p'] 
r(2) 



„ +sgn(7G)Af^ ^ „ 

-p'] 6 V I'^J q,p,p',[q+p-p'] 



(4.9) 

(4.10) 
(4.11) 



Here [q -I- p — p'] stands for q -I- p — p' reduced to 
the first MBZ by a reciprocal vector G, that is, 
[q + p - p'] = q + p - p' - G. Equation and the 

first term of Eq. (|4.9p arise from the spin- reduction factor 
with the diagram (a). The second and the third terms 
of Eq. (I4.9P arise from the diagram (b), and the fourth 
term arises from the diagram (c) . No contribution comes 
out from the diagram (d) . A systematic study of the 1/S 
expansion for the dynamical structure factor has been 
carried out in Ref. 45. 

Figure [5] shows the one-magnon intensity yi{uji]q) 



within the LSW contribution, as a function of q along 
symmetry lines, where uJi is set to give rise to the peak 
in the absorption spectra. Since it is proportional to 
the dynamical structure factor, its q-dependence is well 
known. The one-magnon intensity vanishes with q — >■ 0, 
while diverges with q — (tt, tt). The difference from the 
result of the UCL approximation is only its magnitude; 
the difference becomes larger with decreasing values of 

r. 

Figure [S] shows the one-magnon intensity evaluated up 
to the first and the second order of 1/S [Eqs. (|4.8I) and 
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(|4.9p ]. The intensities are found strongly reduced from 
the zeroth-order values. The l/S expansion has aheady 
been carried out to the transverse component of the dy- 
namical structure factor (see Fig. 5 in Ref. [ish . 



Now we discuss the continuum spectra of three- 
magnon excitations. Analyzing carefully the diagrams 
in Fig. m we finally obtain from the second-order contri- 
bution with respect to 1/5' as 



p.p' 



2'k6{u} - Ep - Cp' - Cq+p-p' 



X -f^ ^2 n2 p2 

^ 2 'I P P 'i+P'P' 



B 



■(4) ss:n('yn)B^°' 
q p.p'Jq-l-p-p'] I ° ^ q:P,p':[q+p-p 



m'-^^ „ +sgn(7G)M*^^ ^ 

q,p,p ,[q+p-p ] ° ^ '^^ q.p.p ,[q+p-p 



(6) 



■^q+p-p 1 



Cq + Sp + Cp' 



(4.12) 



r 



0.2 



0.1 



(a) r/(2J)=1 .2 



• LSW 

1/(2S) 
1/(2S)' 



(b) r/(2J)=0.5 




(0,0) 



(1,0) 



(ii/2,ii/2) (0,0) (0,0) 



(1.0) 



{ii/2,ii/2) (0,0) 



FIG. 6: (Color online) One-magnon intensity 

yi{iJi; ci) /[N{2S)] up to the first (red broken line) and 
the second (green dot-dashed line) order of l/S, as a function 
of q along symmetry lines. The tOi is set to give rise to the 
peak in the absorption spectra. The integrated intensity of 
three-magnon continuum also shown by blue 

thin line. 



symmetry lines for q with uji corresponding to the peak 
in the absorption spectra. It is found that its intensity is 
about 20%-30% of the one-magnon intensity. 

The magnon energy eq is corrected as Cq = (1 + 
A/2S')eq within the first-order in l/S. The second-order 
correction in 1/5*, which is known to be rather small 
(see Ref. [13 for the details). We replace S{uj — Cq) by 
S{uj — €q) in Eq. (|4.6p and S{uj — ep ~ ep> — e^+p^p') 
by (5(w - gp - ep' - e^+p.p/) in Eq. (|4.12p . respec- 
tively. Figure [7] shows the spin-flip correlation function 
y(^) (w^; q, aj)/A^(25') as a function of lj along symme- 
try lines for q. The spectra consist of the (5-function 
peak shown by vertical lines and the continuum of three- 
magnon excitations. 



B. Spin-conserving excitation spectra 



Note that Eq. (I4.12[) is. except for the prefactor, equal to 
Eq. (6.19) in Ref. l45l as it should be. The total intensity 
of three-magnon excitations is given by 



/(3)(w,;q) = y" ?/3(w,;q,cj)^. 



(4.13) 



The two-magnon operator Z^^'>{(jJi;q) defined by 
Eq. (|3.33p is rewritten in terms of magnon operators by 
using Eqs. (jAip - (jA10[) . The result is summarized as 



Z(2)(c.,;q) = (25)^iV(c.,;q,k)at^^^j/3L, 



Figures [6] (a) and (b) show calculated /'"^•'(wi; q) along with 



(4.14) 



N{ujf, q, k) = ^(w,,)|(l + 7q)(^[q+k]"^k + sgn(7G)4m[q+k]) 

+ ^[q+k]4(7k + sgn(7G)7[q+k]) + nT-[q+k]nT-k(7[q+k] + sgn(7G)7k)}, 

I 



(4.15) 



where [q + k] stands for the q -I- k reduced in the first q + k — G. In deriving Eq. (|4.14p . the non-linear terms 
MBZ by a reciprocal lattice vector G, that is, [q -f k] = 
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FIG. 7: Spin-flip correlation function Y''^\uJi;q,uj)/[N{2S)] 
as a function of energy loss uj along symmetry lines for q. 
Vertical solid lines represent the (5-function peak of the one- 
magnon excitation. On panel (b), only one vertical line is 
shown, since the peak positions and the weights of the one- 
magnon excitations are nearly the same for (tt, 0) — (7r/2, 7r/2). 



in the expansion of with magnon operators have been 
neglected, since the consistent analysis within the 1/5 ex- 
pansion is quite complicated. Such terms may cause the 
intensity transfer to the four-magnon excitations. This 
type of correlation function has been studied for the mag- 
netic excitations in the i^T-edge RIXSi ^^i^^ 

Since two-magnons are excited closely to each other 
around the core-hole site, the magnon-magnon interac- 
tion would be important. Introducing the two-magnon 
Green function, 



F(q,c^;k,k') = 



*dt(T[/3_k(t)a 



[q+k; 



q-|-k']^-k'])' 

(4.16) 




(0,0) 



C/3 
— ' 
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FIG. 8: Spin-conserving correlation function 

Y''^''{uji;q,Lj)/[N{2S)^] as a function of energy loss ui 



along symmetry lines for q. 



we take account of scattering of two magnons through 
the term 5^234 [^1- (!A16I) ]. Applying the addition theo- 
rem of trigonometric functions to factors such as 72-4 in 

^1234' could transform -6^234 i^^o a separable form. 
Thereby we obtain the t-matrix by summing up the lad- 
der diagrams in a closed form. We have already explained 
this procedure in the analysis of the two-magnon spectra 
in the K-edge RIXS. See Ref. .25 for the details. Once we 
obtain the Green function, by the fluctuation-dissipation 
theorem, we obtain the correlation function as 

y(2)(^.q,^) ^ (25)2^^7V*(^,;q,k)iV(c.,;q,k') 



(-2)Imi;^(q,c^;k,k'). 



(4.17) 



Figure [5] shows F^^^ (w^; q, w) as a function of uj along 
symmetry lines for q with uji corresponding to the peak 
in the absorption spectra. 

The frequency-integrated correlation function is not 
changed by the presence of the magnon-magnon inter- 
action. Neglecting the interaction in the calculation of 
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(0,0) (7t,0) (7c/2,7t/2) (0,0) 

FIG. 9; Frequency-integrated correlation function 
y'-'^^ (a;;; q)/[A'^(2S')^] as a function of q along symmetry 
lines. The u)i is set to give rise to the peak in the absorption 
spectra. 

F(q, oj; k, k') in Eq. (|4.17l) . we simply obtain 

= (25)2^|iV(^,;q,k)|2. (4.18) 

k 

I 



Figure [9] shows y(2)(^.q)/jv(25')2 as a function of q 
along symmetry lines with uji corresponding to the peak 
in the absorption spectra. The values are about one order 
of mag nitude smaller than those of ^(^^(w,;; q)/[iV(25)]. 



V. SPECTRAL SHAPE IN COMPARISON 
WITH EXPERIMENTS 



In this section, we analyze specifically the Cu L3 edge 
spectra in Sr2Cu02Cl2. According to the experimental 
setup shown in Fig. 1(a) of Ref. 1351, the photon comes 
onto the a-h (x-y) plane, and is scattered with the angle 
130 degrees to the incident direction. The a-c (x-z) plane 
is set as the scattering plane. The polarization vector of 
the incident photon is then expressed as at = (0,-1, 0) 
for the a polarization and oti = (x^, 0, Q) for the tt po- 
larization. Similarly, the polarization of the scattered 
photon is expressed as ct/ = (0,-1,0) for the a' po- 
larization and ckf = (x^,0,(T) for the tt' polarization. 
Thereby we have from Table |lll 



pW(3/2,a^a') = 0, F^^) (3/2, a-^, a') = -, {a ^ a'), 

15 

FW(3/2,a^a^) = -^Xj, P(°)(3/2, a^, a^) = 0, {a ^ n'), 
FW(3/2,a/,aO = 0, (3/2, a') = ^xjx^, (^r -> tt'), 

FW(3/2,a^a^) = ^xl , (3/2, a^, = 0, {tt ^ a'). 



(5.1) 
(5.2) 
(5.3) 
(5.4) 



The polarization is separated with the incident photon, but not separated with the scattered photon in the experiment. 
In this situation, the RIXS spectra, which depend on the polarization of the incident photon, are expressed from 
Eqs. (I5II])-(IM1) as 



AuJiLOf 



yW(a.,;q,c.) + (^)V(2)(a;,;q,c.) 



(f ) V(i) (.;,; q, ^) + ( A) 2 (^^J^^^ ' y (2) (^^ ; q, ^) 



(cr - poL), 

(tt - poL), 



(5.5) 



where q is regarded as the component (of the trans- 
ferred momentum) projected onto the a-b plane. For the 
spin-flip excitations, the intensity in the tt polarization 
is larger (or smaller) by a factor (xf/x/)^ than that in 



the a polarization. For the spin-conserving excitations, 
the intensity in the tt polarization is smaller by a factor 
iXfXi)'^ than that in the a polarization. Notice that, in 
both cases, the intensity ratio between in the tt polariza- 
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FIG. 10: (Color online) RIXS spectra as a function of energy 
loss u for the momentum transfer projected onto the a — b 
plane: panel (a) for q = (0.737r,0), and panel (b) for q = 
(— 0.927r,0). The incident photon energy u)i is set to give 
rise to the peak in the absorption spectra. J = 130 meV, 
and the calculated spectra are convoluted with the Lorentzian 
function with the half width of half maximum 78 meV. The 
experimental data for Sr2Cu02Cl2 are taken from Ref.lH; the 
peak height is adjusted to coincide with the calculated one in 
the TT polarization for q = (O.TStt, 0), while it is adjusted with 
the calculated one in the a polarization for q = (— 0.927r,0). 
The thick solid (black) lines represent the total intensities. 
The thick dotted (red), thin solid (green) and thin dotted 
(blue) lines correspond to the one, two, and three magnon 
contributions, respectively. 



tion and in the a polarization is completely determined 
from the scattering geometry. That is, the ratio is inde- 
pendent of the values of the correlation functions. 

Figure [TU] shows the calculated result with Wi giving 
rise to the peak in the absorption spectra. For compar- 
ison with the experiment, '^^ the exchange coupling con- 
stant J is assumed to be 130 meV, which is nearly the 
same as the one estimated by the inelastic neutron scat- 
tering experiment in La2Cu04.^^ The calculated spectra 
are convoluted with the Lorentzian function with the half 
width of half maximum 78 meV in accordance with the 
experimental resolution. Panel (a) shows the spectra for 
q — (0.737r, 0). According to the experimental setup, we 
have xl — —0.95 and X/ = 0.38. The experimental data 
are drawn such that the peak height coincides with the 
calculated peak in the tt polarization. In the tt polariza- 
tion, the spin-conserving contribution is suppressed by 
the relative weight (2xj)^ = 0.58, and the spectra are 
dominated by the spin-flip contribution. The spin-flip 
contribution includes the three-magnon continuum with 
the energy higher than the single-magnon peak. This 
makes the spectral shape asymmetric in agreement with 
the experimental data in the tt polarization. In the a po- 
larization, the spin-flip contribution becomes smaller by 
a factor (x/Zx^)^ = 0-16 than that in the tt polarization, 
while the spin-conserving contribution becomes larger by 
a factor ^/{x^fXlY = 7.7 than that in the tt polarization, 
and thereby the two-magnon intensity becomes larger in 



comparison with the one-magnon intensity. In the exper- 
iment, the spectra consist of two peaks;^ the low-energy 
peak is considered to come from the one-magnon excita- 
tion, which intensity is estimated to be smaller by a factor 
0.32 than the one-magnon intensity in the tt polarization. 
This value is about twice the theoretical value, and the 
reason for this discrepancy is not known, since the value 
is determined by the geometry. In the calculated spec- 
tra, the one-magnon and two-magnon peaks are closely 
located, forming a single peak. However, the spectral 
shape is asymmetric with a broad width in agreement 
with the experiment, and the total intensity is also in 
good agreement with the experiment. 

Panel (b) shows the spectra for q = (— 0.927r,0). Ac- 
cording to the experimental setup, we have xf — —0.11, 
X} = 0.85. Since {xVx}f = 0.016 and (x^X/)^ = 
0.0087, the intensity in the tt polarization is two order 
of magnitude smaller than that in the a polarization. 
In the fj polarization, the two-magnon intensity has a 
considerable weight, resulting in a shoulder in the high 
energy side. The overall shape with a considerable width 
agrees well with the experiment. 



VI. CONCLUDING REMARKS 

We have studied the magnetic excitations in the L-edge 
RIXS in undoped cuprates. We have analyzed in detail 
the second-order dipole allowed process with paying at- 
tention to the strong perturbation through the interme- 
diate state, in which there is no spin degree of freedom 
at the core-hole site. In this situation, it is not logically 
appropriate to make a perturbation calculation with the 
terms involving the spin degree of freedom lost in the 
intermediate state. Within the approximation that the 
perturbation due to the intermediate state is not extend- 
ing to neighboring sites, we derive the spin-flip final state 
expressed as ck/^ x oti±_ ■ So|.g) in the scattering chan- 
nel with changing the polarization, which leads to the 
RIXS spectra expressed as the dynamical structure fac- 
tor of the transverse spin component. In the scattering 
channel without changing the polarization, we have as- 
sumed a spherical form of the spin-conserving final state, 
X • So Iff), which leads to the RIXS spectra expressed as 
the 'exchange'-type multi-spin correlation function. We 
have numerically evaluated the transition amplitudes for 
both the spin-flip and the spin-conserving final states on 
a finite-size cluster centered at the core-hole site. 

Since no core hole exists in the final state, the spin ex- 
citations could move around the crystal. We have treated 
the itinerant spin excitations by means of the 1/S- 
expansion method, which is known to work for treating 
the quantum fluctuation in the two-dimensional Heisen- 
berg antiferromagneti^i^ For the spin-flip excitations, 
having expanded the spin-flip operators up to the sec- 
ond order of l/S, we have obtained the three-magnon 
excitations in addition to the one-magnon excitations. 
This gives rise to a considerable reduction of the one- 
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magnon intensity as well as the intensity transfer to the 
three-magnon continuum. For the spin-conserving exci- 
tations, we have taken into account the interaction be- 
tween magnons by summing up the ladder diagrams. 

We have analyzed the Cu La-edge spectra in 
Sr2Cu02Cl2 on the basis of these results. The two- and 
three-magnon excitations give rise to substantial intensi- 
ties in the high energy side of the one-magnon peak as 
a function of energy loss, in good agreement with the 
experiment 1^ We hope that the similar analyses are ap- 
plied to the RIXS spectra in other materials and clarify 
the nature of magnetic excitations in future. 
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Appendix A: 1/S expansion 



with ne = ajai and bjbj. Indices i and j refer to sites 
on the up and down sublattices, respectively. Using 
Eqs. (|Al[) - ()A4|) . i/inag may be expanded in powers ofl/S, 

i?„.ag = -^JS^Nz + + + + . . . , (A6) 

where A'' and z are the number of lattice sites and that 
of nearest neighbor sites, respectively. The leading term 

iJmag is expressed as 

<i.3> 

The Fourier transforms of the boson operators are intro- 
duced within the first magnetic Brillouin zone (MBZ), 

/ 2 \ 

ai = i—j E akexp(ik • ri), (A8) 
/ 2 \ i/2 

= [nJ E^kcxp(*k.r,). (A9) 



We summarize briefly the 1/S expansion method in 
the Heisenberg antiferromagnet. For details, see Refs.lsol 
andl45l. The x, y, z axes below are interpreted as the x', 

y', z' 



a. Hamiltonian 

Assuming two sublattices in the antiferromagnetic 
ground state, we express spin operators by boson op- 
erators as 



St = (5r)t = V2^/,(5)a, 



5| = -S 



(Al) 
(A2) 
(A3) 
(A4) 



S+ = (57)t = V2^&]/,(5), 
where and bj are boson annihilation operators, and 



Then, we introduce a Bogoliubov transformation. 



= ^ttk + "^k;5-k, &-k = Tikaj^ + 4^-k, (AlO) 



with 



4 = 



Ek 



l + Ek 



2ek 



1/2 



mk 



Ck 



7k 



2ek 

ik-5 



1/2 



-x^fin) 

(A12) 



where 5 connects the origin with the nearest neighbor 
sites. By this transformation, Hmlg is diagonalized as 



= JSzJ^i^k - 1) + J^z^ ek(4ak + /S^/Sk). 

k k 

(A13) 

Similarly, H^lg is expressed as 



16 



JS 



25* 

—JSz 
2SN 



^(5g(1 + 2-3-4KiW4 



1234 



I 



(6) 
234 



H.c. 



(A14) 



Here the first term in Eq. (|A14p is known as the Oguchi 
corrections^ which coefficient A is given by 



where (■ • • ) denotes the expectation value over the 
ground state, and T is the time-ordering operator. Mea- 
suring energies in units of JSz, the unperturbed propa- 

(0) 



(A15) gators corresponding to H^nlg are given by 



For the square lattice, A = 0.1579. The second term rep- 
resents the interaction between magnons, where ki, k2, 
ks, • • • are abbreviated as 1, 2, 3, • • • , and the Kronecker 
delta 5g(1 + 2 — 3 — 4) indicates the conservation of mo- 
menta within a reciprocal lattice vector G. The vertex 
functions i?*^*)'s in a symmetric parametrization are given 
in Ref. [l^ and [30[ For example, 5(2341 which describes 
the scattering of two magnons, is given by 

(3) 

■^1234 = 72 -4 + 7i-3a;iX2X3a;4 -I- 7i_4a;ia;2 + 72-3a;3a;4 
1 , 

- 2 ('^22^4 + I1X1X2X4 + 72-3-42^3 + 7i-3-4a;ia:2a;3 

+ 742:2 + 732:22:32:4 + 74-2-12:1 + 73-2-12:1X30:4). 

(A16) 

(2) 

The second-order term -ffmag is composed of products 
of six boson operators. Writing it in a normal product 
form with respect to spin-wave operators, we have 



mag 



JSz 



(23) 



2 



+C2(k)(at/3lk + /?_kak) 



(A17) 



Neglected terms are unnecessary for calculating correc- 
tions up to the second order. The explicit forms of Ci (k) 
and C2(k), are given by Eqs. (2.22) and (2.23) in Ref.|M- 



GL(k,c.) 



[uj-ek + iS]-\ (A22) 
Gl^{k,u;)^0, (A23) 
[-L0 - ek + iS]-\ (A24) 



The self-energy is defined by a matrix Dyson's equatio: 



G^AK^) = G%{k,uj) (A25) 
+ EG°^'(k,^)5^MV'(k,c.)G,.,(k,L.). 



It is expanded in powers of 1/(25'), 



I]^.(k,c.) = ^sW(k,c.) + ^E|fj(k,c.) + ■■■. 

(A26) 

The first-order terms are obtained from Hi: 



the Green function 



EL'2(k,w) = 45(k,c^) - Aek, 

Ga/5(k,i) - -i(r(ak(t)/3-k(0))), (A19) 



We introduce the Green functions for spin- waves, \^(i)/i ' '\ -Jl^) ^^ ' \ n ' (A27) 

G„„(k,t) = -i(r(ak(t)4(0))), (A18) 



EW(k,u;) = Eg(k,w) = 0. 

)))), (A18) 
;0))), (A19) 
G^„(k,i) - -i(r(/3lk(i)4(0))), (A20) 



G^/j(k,i) ^ -i{T{l3l^{t)l3^i^{0))), (A21) order perturbation: 



The second-order term S/ji;(k, w) is given by the second- 
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(4) 

k,p,q,[k+p-q] 



(6) 

k,p,q,[k+p-q] 



CJ - ep - eq - ek+p-q + i(5 UJ + ep + €q + Ck+p-q - iS 



(A28) 



ES(k,c.) = C2(k) 



y 244^2^2 sgn(7G)i?['^ rw+ 2(ep + 6q + ek+p-q) 

k p q k+p-q 6 W<^; k,p,q.[k+p-q] k,p,q,[k+p-q] 2 _ I ^ I , ^2 



- (fp + fiq + ek+p-q)2 + iS 
(A29) 



where (5^0, and [k + p — q] stands for k + p — q reduced to the 1st MBZ by a reciprocal vector G, that is, 
[k + p — p'] = k + p — p' — G. In deriving Eqs. (|A28p and (IA29p . we have used the relations, 



-^[k+p-q].q.p.k — ^g^(TG)-Sk.p.q.[k+p-q]) 
^q,[k+p~q],k,p = Sgn(7G)Bk,p,q,[k+p-q]- 



(A30) 
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